Some exact properties of the gluon propagator 



Daniel Zwanziger 1 
'Physics Department, New York University, New York, NY 10003, USA 

Recent numerical studies of the gluon propagator in the minimal Landau and Coulomb gauges 
in space-time dimension 2, 3, and 4 pose a challenge to the Gribov confinement scenario. In these 
gauges all configurations are transverse, d ■ A = 0, and lie inside the Gribov region il, where the 
Faddeev-Popov operator, M(A) = — d^D^A), is positive, that is, (V>, M{A)ip) > for all ip. 

We prove, without approximation, that for these gauges, the continuum gluon propagator D(k) in 
SU(N) gauge theory satisfies the bound i^- / d d k D ^ < N. This holds for Landau gauge, in 
which case d is the dimension of space-time, and for Coulomb gauge, in which case d is the dimension 
of ordinary space and D(k) is the instantaneous spatial gluon propagator. This bound implies that 
limfc^o k d ~ 2 D(k) — 0, where D(k) is the gluon propagator at momentum k, and consequently 
-D(O) = in Landau gauge in space-time d = 2, and in Coulomb gauge in space dimension d — 2, 
but -D(O) may be finite in higher dimension. These results are compatible with numerical studies of 
the Landau-and Coulomb-gauge propagator. 

In 4-dimensional space-time a regularization is required, and we also prove an analogous bound 

on the lattice gluon propagator, ,,} . . f n d d k — ^rSS ^r-rn D ^C^ < N. Here we have taken the 

o x- f o i d{2TT) a J—-K 4 2^ A sin^(fc A /2) — 

infinite-volume limit of lattice gauge theory at fixed lattice spacing, and the lattice momentum 
componant fc M is a continuous angle —it < k^ < tt. Unexpectedly, this implies a bound on the 
high-momentum behavior of the continuum propagator in minimum Landau and Coulomb gauge 
in 4 space-time dimensions which, moreover, is compatible with the perturbative renormalization 
group when the theory is asymptotically free. 



PACS numbers: 12.38.-t, 12.38.Aw, 12.38.Lg, 12.38.Gc 

I. INTRODUCTION 

The successes of perturbative calculations at high en- 
ergy and of numerical studies in lattice gauge theory pro- 
vide strong evidence that the interactions of quarks and 
gluons are correctly described by the non-Abelian gauge 
theory known as QCD. However we lack a satisfactory 
understanding of the mechanism, by which quarks and 
gluons are confined, comparable to that provided by the 
Higgs model of electro-weak interactions. There are sev- 
eral suggestive scenarios which involve the dual Mciss- 
ner effect with condensation of magnetic monopoles, the 
maximal Abclian gauge, the maximal center gauge or the 
light-cone gauge. 

There is also a scenario in Landau and Coulomb gauges 
that originated with Gribov [l| that is based on the in- 
sight that there exist Gribov copies — that is to say 
gauge-equivalent configurations that satisfy the gauge 
condition — and moreover that the dynamics is strongly 
affected if one cuts off the integral over configurations A 
to avoid over-counting these copies. According to this 
scenario, the cut-off is nearby in infrared directions (in 
A-spacc), which suppresses the gluon propagator D(k) 
at small k, so that would-be massless gluons exit the 
physical spectrum and are said to be confined. However 
recent numerical studies, which we shall review shortly, 
have revealed that the behavior of the gluon propagator 
is more complicated than expected, and present a chal- 
lenge to this scenario. In the present article we present 
exact bounds on the gluon propagator, which result from 



the cut-off in A-space, that are consistent with and clar- 
ify the results found in numerical studies of the gluon 
propagator. 

There had been various early conjectures about the 
behavior of the gluon propagator D(k). Gribov in par- 
ticular obtained by an approximate calculation in Lan- 
dau gauge [l[ D(k) = k x +1 in all space-time dimensions. 
This has the notable property that D(0) = 0, in striking 
contrast to the tree- level gluon propagator, D(k) = -p-, 
which has a pole at k = 0. This same propagator is 
also the zeroth-order gluon propagator in a perturbative 
expansion based on a local, renormalizable action that 
includes a cut-off at the Gribov horizon However 
according to numerical studies in Landau gauge it ap- 
pears that -D(O) does not vanish in dimension 2+1 and 
3+1. To address this problem, a dynamically refined ac- 
tion for dimension 2+1 and 3+1 has been proposed and 
studied 0-Q). 

The gluon propagator has been much studied by 
Dyson-Schwingcr equations and related methods. The 
general consensus at present is that there are two types 
of solution, a scaling solution, with D(Q) = 0, and a de- 
coupling solution, with D(Q) > 0, as discussed in ||, 
where further references may be found. 

The gluon propagator in Landau gauge has also been 
the subject of numerical study in lattice gauge theory. 
From recent studies, a somewhat puzzling picture has 
emerged. It appears that in dimension 1+1 the gluon 
propagator in Landau gauge does in fact vanish at k = 0, 
D(0) = 0, in accordance with Gribov's result [9j-|Li|. In 
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dimension 2+1 it was found that the gluon propagator 
has a turn-over below which D{k) decreases with decreas- 
ing fc, and there is no explanation for this non-intuitive 
behavior besides the proximity of the Gribov horizon in 
infrared directions. However from studies on huge lat- 
tices, it appears that in dimension 2+1, D(k) approaches 
a, finite value at k = [HI EH, [H| . In dimension 3+1 there 
may be a kind of shoulder at low momentum and, as in 
dimension 2+1, it appears that D(0) is finite [II [ll, GJ- 

In Coulomb gauge the numerical results for the equal- 
time space-space propagator D(k) in space-time dimen- 
sion d + 1 are qualitatively similar to those in Landau 
gauge in space-time dimension d, and infrared suppres- 
sion is more pronounced for d — 2 than for d = 3, with a 
more distinct decrease with |k| at small k for d = 2 fl~7i — 
|2~H . 1 Most of these numerical studies of Coulomb gauge 
are for SU(2) gauge theory in dimension 3+1, with re- 
sults for dimension 2+1 in |19|. [2ll and for SU(3) gauge 
theory in dimension 3+1 in |20(. However from numer- 
ical studies there is as yet no definite conclusion as to 
whether in Coulomb gauge -D(O) vanishes in dimension 
2+1 and 3+1. 

To summarize, we are faced with the puzzle that nu- 
merical studies in Landau gauge indicate that D(Q) = 
in dimension 1+1, which accords with the Gribov sce- 
nario, but -D(O) is positive, D(0) > 0, in dimension 2+1 
and 3+1. The puzzle deepens in view of an argument 
that led to the conclusion that D(Q) = in any number 
of dimensions (22[. This argument involves the free en- 
ergy, W(J) = \n(e( J ' A )} : where J is an external source. 
The input for this argument is (i) a bound on W(J), and 
(ii) the hypothesis that W(J) is analytic in J. The bound 
on W{J) appears unassailable, and so, if the numerical 
results are accepted, it must be that W( J) is non-analytic 
in J, which is the signal for a change of phase. The non- 
analyticity of W(J) will be discussed elsewhere J23|. 

In the present article we take up another bound, which 
is found in Appendix B of [22[ , that does not involve the 
free energy W(J). It is an ellipsoidal bound, satisfied by 
all configurations A in the Gribov region 57. From this 
bound on configurations A, we will obtain a bound on 
the gluon propagator D(k) by taking expectation values. 

In sect. II A we recall some known results, in sect. IIB 
we explain the simple idea on which the derivation of the 
bound is based, and in sect. IIC we note that the bound 
applies to other gauge bosons. In sect. Ill we present 
exact bounds on the gluon propagator that hold in con- 
tinuum and in lattice gauge theory, and we also present 
the renormalized continuum bound that holds in dimen- 
sion 3+1. In sect. IV we discuss the implications of these 
bounds for the Landau- and Coulomb-gauge propagator 



in dimension 1+1, 2+1 and 3+1. In addition to infrared 
bounds that are stronger in lower space-time dimension, 
we shall find, unexpectedly, that the cut-off at the Gribov 
horizon implies a bound on the /izg/i-momentum behav- 
ior of the Landau- and Colulomb-gauge propagator in 
dimension 3+1 given in eqs. (|35|) and (|50|) respectively. 
Some concluding remarks may be found in sect. V. In 
Appendix A we derive an ellipsoidal bound on contin- 
uum configurations that lie inside the Gribov horizon. 
In Appendix B we convert this into a bound on the con- 
tinuum gluon propagator. 2 In Appendix C we exhibit a 
simpler ellipsoidal bound on continuum configurations in 
the infinite- volume limit. In Appendix D we derive the 
bound on the lattice gluon propagator. 

II. SET-UP AND BASIC IDEA 

A. Elementary properties 

We deal with Euclidean QCD in its continuum and lat- 
tice formulations. Numerical gauge fixing on the lattice 
is done by gauge transforming to a local minimum of a 
lattice analog (specified in Appendix D) of the continuum 
minimizing functional |24l426l | 

F A {g) = I M| 2 

= J d d x \ 9 A(x)\ 2 . (1) 

This is the Hilbert square norm of the configuration 
9 ' An, = g~ 1 A^g + g~ 1 d ll g, which is the gauge-transform 
of the gauge field A^x) = t b A b ^(x) by the local gauge 
transformation g{x) S SU(N), where fi — l,...d. Here 
the t a are an anti-hermitian basis, (t a ) ' = —t a , of the Lie 
algebra of the SU(N) group, [t a ,t b ] = f abc t c , normalized 
to tr(t a t b ) = -S ab /2, where a = 1, ... N 2 - 1. If d is 
taken to be the dimension of Euclidean space-time, then 
this gauge fixing produces a gauge in the class of minimal 
Landau gauges, whereas if d is the dimension of ordinary 
space, then the gauge is in the class of minimal Coulomb 
gauges. (For the Coulomb gauge this minimization is 
done at every Euclidean time t.) The minimization pro- 
duces a local minimum of the minimizing functional. Any 
local minimum will do. In principle it could be the abso- 
lute minimum, but this is not necessary for our purposes, 
nor is it achievable in practice numerically. 

At a local minimum (i) the functional Fa(9) is sta- 
tionary, and (ii) the matrix of its second derivatives is 
positive. Property (i) gives the transversality condition, 

MS = (2) 



1 Reference [2l| also presents results for gauges that interpolate 
between Coulomb and Landau gauges, with gauge condition 
XdoAo + diAi = 0. 



2 The same continuum bound is also obtained in cq. (1 1 5 1) as a limit 
of the lattice bound derived in Appendix D, but we have provided 
an independent derivation of the continuum bound because it is 
simpler. 
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characteristic of the Landau gauge. Property (ii) is the 
positivity of the Faddccv-Popov operator, 

a>,M(A)i>) = (d^,d^) - (rr bc ,A%r) > o, (3) 

for any wave function ip a (x) . These two properties define 
the (first) Gribov region f2, and gauge fixing by this min- 
imization produces configurations A that all lie inside fi. 

B. Basic idea of bound 

It is very easy to establish bounds on configurations 
A that are in the Gribov region ft. Take any trial wave 
function ij)(A) that may depend on A. Then, from it 
follows that every A in fi satisfies the bound, 

(r(A)f abc ,A%ip c (A)) < (d^(A),d^(A)). (4) 

For an appropriately chosen trial wave function ip(A), an 
ellipsoidal bound on A of the form, 

£ c &X>m<i. ( 5 ) 

k 

is obtained, as shown in Appendix A. Here a\ is the 
component of A' {x) in the fourier expansion, 

A b ll (x) = J24,^ ik - x , (6) 
k 

on a finite periodic Euclidean volume, V = L d , where 
kfj, = 2irn fl /L, and runs over all integers. Such a 
bound for a finite lattice was established in Appendix B 
of (2^ | , and a stronger ellipsoidal bound is derived in the 
present article for continuum and lattice gauge fields in 
Appendices A and D respectively. Upon taking expecta- 
tion values, we obtain the bound on the gluon propagator 

v~ x Y, c£(*)JV(*) < 1, (7) 

k 

where we have used (a£* a£ „) = V^ 1 S bc D lll ,(k). In di- 
mension 3+1, the continuum theory must be regularized, 
and in Appendix D, a bound on the lattice gluon propa- 
gator is derived from the positivity of the lattice Faddeev- 
Popov operator. The limit in which the ultraviolet regu- 
lator is removed, A — > oo, is discussed in sect. IV. 

C. Other gauge bosons 

The only input to the bounds obtained here is the re- 
striction of the functional integral to the interior of the 
Gribov region. For this reason, the bound is the same 
whether or not the gluons are coupled to quarks or not, 
although the bound becomes inconsistent in dimension 
3+1 if the theory is not asymptotically free. The bounds 



obtained here also apply to the propagator of other gauge 
bosons that belong to an SU(N) gauge group, includ- 
ing those with Higgs coupling. In the present article we 
are concerned with QCD gauge theory only. However 
it should be noted that the Landau gauge is a special 
case of the gauge, with £ = 0, that is used when the 
gauge field is coupled to a Higgs boson. This gauge may 
be given a non-perturbative meaning by the minimizing 
gauge fixing described above. This is straightforward in 
dimension 1+1 and 2+1, and our results hold in these 
dimensions. In dimension 3+1, a lattice regularization 
of ultraviolet divergences would be required to give the 
theory a non-perturbative meaning, but we have not con- 
sidered other gauge bosons in dimension 3 + 1. 



III. BOUNDS ON GLUON PROPAGATOR 
A. Bound on continuum gluon propagator 

The continuum propagator is defined by 

(Al(x)Ai(O)) = V- 1 £ S bd D^(k)e ik - x , (8) 
k 

where the (hyper)cubic periodic volume V = L d is suffi- 
ciently large that (hyper) spherical symmetry holds, and 
the propagator has the tensor structure 

D ltv (k) = D(k){s itu -^-) f (9) 

by transversality. Note that a factor of the coupling con- 
stant go has been absorbed into the gauge field A, so 
D = cjqDq, where Dq is the unrenormalized, canonical 
propagator. 

Statement: For gauge fixing to the interior of the Gri- 
bov region f2, as in the minimization procedure described 
above, the gluon propagator D(k) in SU(N) gauge the- 
ory satisfies the bound 

d (2n) d J k 2 ~ K ' 

This holds for Landau gauge, in which case d is the di- 
mension of space-time, and for Coulomb gauge, in which 
case d is the dimension of ordinary space and D(k) is 
the instantaneous spatial gluon propagator. (Because of 
renormalization, dimension 3+1 requires a special dis- 
cussion which will be given shortly.) This is proven in 
Appendices A and B. 

Lest it be thought that this bound is trivial, note that 
for the free propagator in Landau gauge, D(k) = l/k 2 , 
the bound is violated because of an infrared divergence 
of the k integration in dimension 1+1, 2+1, and 3+1, 
and by an ultraviolet divergence of the k integration in 
dimension 3 + 1 . 
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If the angular integration is performed, the bound 
reads 



= d . 1 S *-\ I d k k d ~ 3 D(k) < N, 



d (2ir) d 



where Si = 2tt, S2 = 47r, and S3 = 2tt 2 . Note that J has 
engineering dimension in all dimensions d. 



B. Bound on lattice gluon propagator 

To discuss the case of dimension 3+1, we must regu- 
larize and renormalize. The lattice provides a convenient 
regularization, and a lattice analog of the bound fflO|) 
holds, 



the lattice bound §F2§ goes over to the bound, 

CA 



J 



1 Sd-i 



d (27r) d 



dk k d ~ 6 D{k) 



(15) 



,9o 2 (A)^ 3 (A) 



d-l S d -i 



CA 



d (2 



*-± J dk k d - 3 D r (k) < N. 



This is the renormalized form of the continuum bound 
(fTTI) that holds in dimension 3+1, where d = 4 for the 
Landau gauge and d = 3 for the Coulomb gauge. We will 
find that the limit A — > 00 is independent of C. 



IV. DISCUSSION 



A. Infrared bound in Landau Gauge 



d(2ir) d 



dd/; E M cos 2 fa/2) D^{k) 
4E A sin 2 (fcA/2) 



<N, (12) 



as is shown in Appendix D. (We have used the same 
symbols J.D^... for continuum quantities and their 
lattice analogs.) Here we have taken the lattice 
volume to infinity keeping the lattice spacing finite, 
and the lattice momentum fc M is a continuous angle 
—it < fc M < 7r. The lattice propagator is given by 
D^(k) = T, x (^ b x,^o,u)^ iHx+(efl/2) ~ ie '' m] , where the 
lattice variable A\ is defined in (|D3[) . 



Renormalized form of the continuum bound 



The lattice variable A goes over in the continuum limit 



to 



A -> g (A)A Q 

= g {K)Zl /2 {K)A r , 



(13) 



where go = <?o(A) is the bare coupling constant that de- 
pends on the ultraviolet cut-off A, Aq is the unrenormal- 
ized or canonical continuum gauge field, and A r is the 
renormalized continuum gauge field. Consequently the 
lattice gluon propagator is related to the unrenormalized 
and renormalized gluon propagators by 



D -> g 2 (A)D 

= g 2 (A)Z 3 (A)D r , 



(14) 



where D r {k) is the finite, renormalized, continuum prop- 
agator. The lattice momentum goes over to the contin- 
uum momentum by k — > ak, where a is the lattice spac- 
ing, so the lattice integral /_ d d k goes over to the contin- 
uum integral with a cut-off J^ <c ^d d k 7 where A = 1/a, 
and C is a constant of order 1. Consequently, for large A 



We first discuss the Landau-gauge case. The bound 
pip is more stringent in the infrared in lower dimensions 
because of the factor k d ~ 3 in the integrand. Since the 
bound is finite and the integrand is positive, it follows 
that in dimension d = 1 + 1 the propagator D{k) must 
vanish at k = 0, 



lim D(k) 

fc-s-0 



for d = 1 + 1. 



(16) 



However in dimension 2+1, the bound (|11|) is compatible 
with a finite value for D(0). Interestingly that is precisely 
the behavior that has been found in numerical studies of 
the gluon propagator in Landau gauge in dimension 1 + 1 
and 2 + 1, as discussed in the Introduction. D(k) may 
even be singular at k = for d = 2 + 1, provided that 
the strength of the singularity remains less than 1/k, 



lim kD(k) 

fe^O 



for d = 2 + 1. 



(17) 



This condition forbids the existence of gluons of mass 
zero for d = 3. 

Numerical studies also indicate that in dimension d = 
2 + 1 the propagator in Landau gauge is suppressed in 
the infrared, although not as severely as in dimension 
d = 1 + 1, with D(k) decreasing with k as k decreases 
to 0, but approaching a finite value, -D(O) > 0. There is 
no other explanation for this otherwise counter-intuitive 
decrease besides the proximity of the Gribov horizon in 
infrared directions. The bound (fTT|) for d = 2 + 1 does not 
require such a decrease, but only that D(k) not diverge 
as strongly as 1/k. Thus it appears that this bound by 
itself does not fully express the strength of the dynamical 
consequences of the cut-off of the functional integral at 
the Gribov horizon in dimension 2+1. 

According to the lattice bound (fT2 
dau propagator D^ik) in dimension 
a singularity as strong as ^ 



sin 2 (fe M /2) 



, the lattice Lan- 
3+1 cannot have 
at k = (where 



cosfa/2) = 1), 



lim sm 2 (V2) ^D vu {k) = 0. 



(18) 
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Thus the bound on the lattice gluon propagator in Lan- 
dau gauge in 4 Euclidean dimensions does not tolerate a 
l/sin 2 (fc/2) singularity, that is to say, a massless lattice 
gluon, for any finite value of the lattice spacing a or, in 
other words, for any finite value of the cut-off A = 1/a. 
Moreover the finiteness of the renormalized continuum 
bound (|15[) for large but finite A implies 



lim k 2 DJk) = 0, 

fe->0 



for d = 3 + 1. 



(19) 



According to numerical studies discussed in the Intro- 
duction, in dimension 3+1 there appears to be a kind of 
shoulder in D(k) at low k, with a finite value of D(0) > 0. 
This is entirely consistent with the infrared bound ob- 
tained here, but again, as in dimension 2+1, the gluon 
propagator D(k) is apparently more strongly suppressed 
in the infrared than required by the bound we have ob- 
tained. Thus it appears that also in dimension 3+1, the 
bound obtained here does not by itself fully express the 
strength of the dynamical consequences of the cut-off of 
the functional integral at the Gribov horizon. 



B. High-momentum bound in Landau gauge 

The coefficient Qq(A)Zs(A) which appears in the renor- 
malzied continnum bound (|15p is either zero or infinite 
as A — > oo, so it might be thought that if <?q(A)Z3(A) 
is zero, then the bound (fl~5|) is trivially satisfied by any 
finite renormalized propagator D r , and if <7g(A)Z3(A) is 
infinite, then the bound implies that the renormalized 
propagator D r (k) vanishes identically for all fc, in which 
case the theory is inconsistent. However before coming to 
this conclusion, we must also consider the A-dependence 
introduced by the ultraviolet cut-off of the integral at 
k = CA. 

In dimension d = 3 + 1, the bound (11511 reads 



CA 

j ee 9 2( A)Z3 ( A )_jL r dk kDr ( k) < K 



(20) 



We evaluate J using the perturbative renormalization- 
group according to which, asymptotically at large A, 



9o 2 (A) 



1 



2b In A ' 



(21) 



where b is the leading coefficient of the /3-function 
dgo/d\nA = — 6o 3 + 0(g 5 ). It is gauge- independent and 
has the value [231, P- 653, 



1 /UN 2n f \ 
(^\~3 3~J' 



(22) 



where n/ is the number of quarks in the fundamental 
representation. The dependence of Z 3 (A) on A is found 
from the perturbative renormalization group. We have 



l + c^lnA + O^) 



(23) 



where, in Landau gauge for SU(N) gauge theory J27|, p. 
589, 



1 



(4tt) 2 V 3 



/13N 4n f 



(24) 



and we have added the quark contribution. According to 
the renormalization group we have 



tf In Z-i 



d\nA 



= cg r 



0(g 4 r ) 



cgl + 0{gt) 
c 



2b In A 



which is solved by 



Z 3 (A) = z 3 (lnA)P, 



(25) 



(26) 



where 



P = 2b 



137V - 4n 



22N - 4n/ ; 



(27) 



and Z3 is a finite constant of integration. Thus the renor- 
malized bound at large A reads 



3^3 

64tt 2 & 



CA 



ln p_1 A J dk kD r (k) < N. 



(28) 



If p > 1, the coefficient of the integral, ln p_1 A, diverges, 
so the bound, J < N, requires that the renormalized 
propagator, D r {k) 1 vanish for all fc, in which case the 
theory is inconsistent. 

To see if this happens, we note from (|27|) that 



p< 1, 



(29) 



provided that the denominator is positive, that is, pro- 
vided that lliV > 2nf. This is the restriction on the 
number of quarks for the theory to be asymptotically 
free, b > 0. Thus the inconsistency is avoided provided 
that the theory is asymptotically free, as we now assume. 

Because we are in the case p < 1, the coefficient of the 
integral in ([28]) vanishes in the limit A — > oo, 



(30) 



lim ln p_1 A = 0. 



There are three possibilities, (i) If the integral in 
is finite, the bound is trivially satisfied — and vacuous. 
This also happens if the integral diverges at large A, but 
too weakly to compensate the vanishing of the coefficient, 
(ii) If the integral diverges sufficiently strongly that J di- 
verges for A —> oo, then the theory is inconsistent, (iii) If 
the integral has a divergence that precisely compensates 
for the vanishing of the coefficient at large A, then a finite 
bound results. 
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To find out which possibility is realized, we first note 
that any finite contribution to the integral is annihilated 
by the coefficient, and we may write the bound as 



CA 



ln p_1 A / dk kD r (k) < N, (31) 



3^3 

64tt 2 6 



where \i is an arbitrary mass. Only the asymptotic form 
of D r (k) at large k concerns us. 

According to the Callan-Symanzik equation, the cor- 
rections to scaling in the gluon propagator are loga- 
rithmic at high momentum, and the renormalized gluon 
propagator has the asymptotic behavior 



D r (k) 



fc 2 ln p k ' 



(32) 



where the same power p appears here as in Z3, and r is 
a finite constant that depends on the normalization con- 
dition. Since any finite contribution to the integral gets 
annihilated by the vanishing of the coefficient, ln p_1 A for 
A — >• 00, we may extend this asymptotic expression down 
to finite k, and the quantity J in (|31|) is thus, asymptot- 
ically at large A, given by 



J = 



3z 3 r 

64tt 2 6 
3z 3 r 



CA 



h/- 1 A 



In" -1 A 



dk/k 

]n 1 - p (CA)-]n 1 - p ii 



647r 2 6 1 — p 

and, with p < 1, we obtain the limit 

lim J = — ^5— 

A->oo 647T^O 1—p 



(33) 



3z 3 r 



1 



(34) 



32tt 2 2b -c 
2N' 

The result is independent of [i and C, and finite. From 
liniA_ s . 00 J < N we obtain the non-trivial bound, 

z 3 r < 2N 2 . (35) 

Four comments: (i) The trivial and inconsistent possi- 
bilities are avoided because the divergence of the integral 
compensates for the vanishing of the coefficient in the 
limit A — > 00. 

(ii) The bound on the high-momentum limit of the 
gluon propagator is rather unexpected because the only 
input to this bound is the restriction of the functional in- 
tegral to the interior of the Gribov region, and it was gen- 
erally believed that this has dynamical consequences only 
in the infrared. Note however that the high-momentum 
bound occurs only in dimension 3+1. 

(iii) The renormalized bound (f5U|) may be expressed in 
terms of unrenormalized quantities 



CA 

32^ / dkk D ^ = I ^ 



(36) 



Since the left-hand side is constructed out of unrenor- 
malized quantities, it is independent of any renormaliza- 
tion scheme. On the other hand it is finite and indepen- 
dent of A at large A. As such, it is a renormalization- 
group invariant. Thus the finite quantity z 3 r for which 
we have just established the bound z%r < 2N 2 , is in fact 
a renormalization-group invariant, although this was not 
apparent from the way Z3 and r were introduced. 

(iv) The number n/ of quark flavors has dropped out 
of the inequality (f3"5"j) . 



C. Infrared bound in Coulomb gauge 

The minimal Coulomb gauge is obtained by minimizing 
the Hilbert square norm of the space components A^(x, t) 
on each time-slice t, 

F A (g,t) = \\ 9 Mt)\\ 2 = f d d x \ 9 M*,t)\\ (37) 



where i = 1, ...d, and the dimension of space-time is d+1. 
Consequently the equal-time propagator 

£>(k)(<y y -fefc i /k 2 )5 ic = J d d x e-**(4(x,t)A|(0,i)), 

(38) 

satisfies in space dimension d the bounds we have 
derived in Landau gauge in space-time dimension d. 
The expectation-value is independent of t by time- 
translation invariance. For orientation purposes we note 
that in zeroth-order perturbation theory the equal-time 
Coulomb-gauge propagator is given by 



D<°>(k) 



dko 1 



2tt k 2 +k 2 2|k|' 
For space dimension d = 2 we obtain from (jlll) , 



lim D(k) = 0, for d = 2, 

k-s-0 



(39) 



(40) 



so the equal-time propagator vanishes at k = 0, D(0) = 
0. This states that a gluon of zero momentum cannot 
be created by applying the field A^(x,t) to the vacuum, 
and thus the would-be physical gluons exit the spectrum. 
From (|T5l) . we obtain, as in (fl~9l) . 



(41) 



lim |k|A-(k) = 0, ford = 3, 

k— >0 



where D r (k) is the renormalized, equal-time, space-space 
Coulomb-gauge propagator. 



D. High- momentum bound in Coulomb gauge 

For space dimension d = 3, we must consider regular- 
ization and renormalization, as in the Landau-gauge. Al- 
though renormalization in Coulomb gauge has not been 
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established to all orders, we suppose that it is renormal- 
izable, and that the perturbativc rcnormalization-group 
holds. 

We proceed exactly as in the Landau gauge case, but 
with different values for the constants. For space dimen- 
sion d = 3, eq. (fT5|) reads 



CA 



J^. 9o 2 (A)Z 3 (A)^ J dk D T (k) 



<N, 



(42) 



where go (A) is the unrenormalized coupling constant, and 
Z 3 (A) is the renormalization constant for the space com- 
l li 

ponents, A^o = Z s A^ T , in Coulomb gauge. (The apace 
and time components of renormalizc differently in 
Coulomb gauge.) 

The renormalization constant Z 3 of the space compo- 
nents Ai of the gauge field is given by (f23|) where, by 
cq. (B.37) of [H], with Z 3 = Z\, the coefficient c has the 
value, 



(4tt) 



1 4rif 



(43) 



for SU(N), and we have added the quark contribution. 
As in the Landau gauge, we have 



Z 3 (A) = z 3 ln p A, 



(44) 



where p — c/2b. Here c, r and z 3 have values appropriate 
to the Coulomb gauge, and we obtain 



CA 



(45) 



where the gauge- independent quantity b is given in (122p . 
As in Landau gauge, the theory is consistent only if p < 1. 
We find 



c _ 3N-2n f 
P ~ 2b ~ UN -2n f '' 



(46) 



and, as in Landau gauge, we have p < 1, provided that 
the denominator is positive, UN > 2nf. This is, again, 
the condition on the number of quarks for the theory 
to be asymptotically free, as we now assume. As before, 
the coefficient of the integral in (|45|) vanishes for A — > 00, 
and any finite contribution to the integral is annihilated 
in this limit. 

According to the Callan-Symanzik equation, the renor- 
malized equal-time propagator has logarithmic correc- 
tions asymptotically at large k, 



£> r (k) 



2|k|ln p |k|' 



(47) 



where p is given in (|46[) . and, as in (|33p. we have 



CA 

= _£i_ hi?- 1 A f dk - I „ , < N. (48) 



6ir 2 b 



2k\n p k 



As in Landau gauge, this gives, asymptotically at large A, 

z 3 r 



J = 



We thus obtain 



12tt 2 6(1 -p) 
6tt 2 (26- c) 
2JV" 



z 3 r < 2N J 



(49) 



(50) 



the same bound as in Landau gauge and, again, the num- 
ber of quark flavors has dropped out. 

As in Landau gauge, the quantity z 3 r is a 
renormalization-group invariant. In Coulomb gauge the 
bound on this quantity governs the high-momentum limit 
of the space-components of the equal-time gluon prop- 
agator -D(k), whereas in Landau gauge it governs the 
high-momentum limit of the Lorentz-invariant propaga- 
tor D(k). 



V. CONCLUDING REMARKS 

We have obtained the continuum and lattice 
bounds (fTTj) and (fl~2j) on the gluon propagator, and the 
bound (|15p on the rcnormalizcd gluon propagator which 
hold in Landau gauge, where d is the dimension of space- 
time, and in Coulomb gauge, where d is the dimension of 
space and D(k) is the instantaneous spatial gluon propa- 
gator. In space-time dimension 2, 3 and 4, these bounds 
imply restrictions on the infrared behavior of the con- 
tinuum gluon propagator in Landau and Coulomb gauge 
that are more severe in lower dimension, and in space- 
time dimension 4 there is, unexpectedly, a restriction 
on the high-momentum behavior of the continuum gluon 
propagator in Landau and Coulomb gauge. 

It would be of interest to test the lattice and contin- 
uum bounds using numerical lattice data for the gluon 
propagator in 2, 3, and 4 space-time dimensions in Lan- 
dau and in Coulomb gauge. It is possible that the bounds 
are not close to being saturated. 
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Appendix A: Ellipsoidal bound on continuum 
configurations 

We shall establish an ellipsoidal bound on continuum 
configurations A that lie inside the Gribov region VI. We 
consider the SU (N) gauge group on a periodic Euclidean 
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volume V = L d . To start, we substitute into inequal- 
ity (|4]) the trial wave function 



ip(A) =ip — a 



Pn 



(M -p 2 ) 



M^Atyo, (Al) 



where M = M +Mi is the Faddecv-Popov operator, with 
M = -d 2 , M? C {A) = -f^A^d^. This wave function 
is inspired by first-order perturbation theory, according 
to which the first-order change in the zcroth-order wave 
function ipo is given by a similar expression, but we shall 
of course obtain an exact bound. The plane-wave state 



ipo = V 1 I 2 exp(zp ■ x) r), 



is an eigenvector of Mo, Moipo = p 2 ipo> rj is an x- 
independent, normalized color-vector, and p^ is an al- 
lowed momentum vector on the periodic Euclidean vol- 
ume V = L d , = 2irn^/L, where is an integer. The 
operator P q is the projector defined by the kernel, 



P q (x, y) = V^ 1 



E 

k ; | k | > | q\ 



Jk-(x-y) 



which projects onto the direct sum of eigenspaces of Ma 
belonging to all eigenvalues k 2 of Mq that are greater 
than or equal to some fixed eigenvalue q 2 . We stipulate 
that the inequalities, 



k 2 >q 2 >p 2 > (2tt/L) 2 



(A4) 



are satisfied so the denominator in (|A1[) is always posi- 
tive, and that k^^p^, and are allowed momentum vec- 
tors on the periodic volume V = L d . The quantities a, 
p^ and q^ are at our disposal, and a will be a varia- 
tional parameter. Note that ipo is independent of A, and 
Mi = Mi(A) is linear in A, so the trial wave function 
%jj = ip(A) has a piece that is independent of A and a 
piece that is linear in A. 

With this wave function we have, by ([!]), for all A e fi, 



(V>, Mtp) = 1(a) = X- 2aY + a 2 Z > 0, 



(A5) 



where 



X 
Y 

Z 



(ipo, (M + Mi)^o) 



(V'o.Mi 



(M -p 2 ) 
(M -p 2 ) 



M^o) 



x(M +Mi) 



(M -p 2 ) 



M^o), (A6) 



and we have simplified Y using P q Moipo = 0. Positivity 
of 1(a) for all a implies 



X > 0; 



Z > 0. 



(A7) 



Moreover /(a) has a minimum at a = Y/Z, from which 
we obtain the bound 



Because A appears only in Mx(A) which is linear in A, 
there is a term in Z which is cubic in A, whereas X and 
Y are at most quadratic in A. We shall bound the cubic 
term by means of the following lemma. 
Lemma: For the SU(N) group the bound 



(w, [M + Mi(A)]w) < N 2 (uj, M uj) 



(A9) 



holds for any A in H and any wave function cu(A). 
This lemma is derived in Appcnix B of [22| . but we 
present the derivation here for completeness. 
Proof: Consider first the SU(2) group. We decompose 



(A2) Mi into the sum of 3 operators, Mi = M x (1) +M 



r(2) 



-M 



(3) 



that are each given by 



M[ b) (A) = S b H b (A) 



(A10) 



(no sum on b), where (S ) a 



b\ — 



ibe 



is an angular mo- 



mentum matrix in the spin-one representation that acts 
on color variables, and H b = iA b is a Hcrmitian opera- 

(A3) 

tor that acts only on space variables. We first bound the 
operator m[ 3) = S 3 H 3 . Let = e%<f>(x), where <j)(x) 
is any function of x, and e± and eg are x-independent, 
normalized eigenvectors of S3, S^e m = me m . We have 

< (o>±, [M + Mi(A)]w±) = (0, M o 0) ± (0,-ff 3 (A)0), 

(All) 

where the inequality holds for all A € SI, and we have 
used the fact that (e±,S 1 e±) = (e±,S 2 e±) = because 
Si and 52 are off-diagonal in the S3 basis. It follows that 
the inequality 



\(^H 3 (A)<j>)\ < (4, Mo® 



(A12) 



holds for any cj>(x) and all A € SI. We now decompose any 
wave function u> according to cu = e+4>+ + eo4>o + e -0- 5 
and we have for all u> and all A <G SI 

\(uM 3 \A) U )\ = \{u,S a H 3 u>)\ 

= |(0 + ,tf 3 </>+)- W>-,# 3 <MI 
< |(0 + ,tf 3 </>+)l + IW>-,# 3 <MI 



< 



,M O 0+) + (</>-, Mot-) 



< (w,Mqw). 



(A13) 



The same inequality holds for and M^ which gives 



|(w,Mi(,4)u;)| < 3(w,M w). 



(A14) 



For the SU(N) group, the proof is identical except that 
there are N 2 — 1 terms in Mi, which gives 



|(w,Mi(A)w)| < (N 2 - l)(w,M w), 



(A15) 



and (|A9|) follows □. 

We apply the lemma to Z which is of the form Z = 



(w, MuS), where u> 



(Mo-p 2 ) 



Mi -0Q. The lemma yields 



Y 2 < XZ. 



(A8) 



Z < N 2 Zo 



(A16) 
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where 



P A/f 

Zo^^M^A) — q 2 Mi(A)^o), (A17) 



(A/0-P 2 ). 
and we obtain the bound. 

Y 2 < N 2 XZ 



(A18) 



The gain here is that Zq is only quadratic in A, whereas 
Z contains a term that is cubic in A. 

We further simplify the bound by comparing Y and 
Zq. We insert a complete set of eigenstates, 

^k,a = V-^e^ea, (A19) 

where e a are a basis of color vectors, and obtain 

1 

k 2 

k . a ; | k | > | q \ 

E 



y= E 



fc 2 -p 2 
(k 2 -p 2 ) 2 



U k ^MM 2 (A20) 
IW>M> Wo)| 2 - (A21) 



From the restriction k 2 > q 2 > p 2 it follows that 



k 2 q 2 
< — — 



k' 2 — p 2 q^ — p^ 



and consequently 



Zq < 



q 2 — p 2 



Y. 



This gives the bound 



q 2 N 2 q 2 
ZqY < ttY" < tjXZq, 



<7 Z — p z 



q z — 



and so 



(A22) 



(A23) 



(A24) 



(A25) 



We next consider 



X = (V> , (M + MOVo) 

= P 2 - ip».f abc V a *a h k=0 ^, 



(A26) 



where we have used the fourier expansion Q , and a b k=0 
is the 0-momentum componant of A h ^{x). We next show 
that 



fabc a* b c 



< 



2n 

T' 



(A27) 



where /i is a fixed Lorentz index. To do so, we 
use (u>, M(A)lu) > for A G ft, where oj = 
V" 1 ! 2 exp(±i2-!Tx Pi /L)r]. We have 

(w.Mw) = (^) 2 T »(j)f k fi„^ c > 0, (A28) 



from which (|A27[) follows. This gives 

X <p 2 + (2tt/L)Y,\p, 

and we obtain the bound on Y, 
N 2 q 2 



Y < 
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P 2 + (2tt/L)J2\p, 



(A29) 



(A30) 



Note that Y is quadratic in A, while the right-hand side 
is independent of A, so this is an ellipsoidal bound on A, 
as advertised. 

We next evaluate Y, eq. (|A20j) . We have 



(V>fc,a,MiV> ) 

This gives 

Y 



-iPuf^V- 1 J d d x Al(x) 



(A31) 



E 

k.a: \ k\>\q\ 



\pj abc a b k _ p ^ c \ 



el 2 



E 

k,a;\k-\-p\>\q\ 

and we have the bound 



k 2 — p 2 



Wf abc 4,X\ 2 

(k + p) 2 — p 2 



r^rW p»pv E 

k; | k-\-p \ > \ q 

N 2 q 2 



< 



q 2 — p 2 



(k+p) 2 -p 2 

P 2 + (2Vi)Ei^ 



(A32) 



(A33) 



For each p and q satisfying q 2 > p 2 > (2ir/L) 2 , and for 
each color vector 77, this is an ellipsoidal bound on the 
fourier componants a k that holds at finite Euclidean 
volume V for all configurations A £ fi. Geometrically 
speaking, the configurations A that satisfy the bound 
()A33l) define an ellipsoid E in configuration space (that 
depends on p and q) . The Gribov region f2 is contained 
in E, and we have the inclusions 



AancE. 



(A34) 



Here A is the fundamental modular region which consists 
of the absolute minimum of the minimizing functional on 
each gauge orbit. 



Appendix B: Bound on continuum gluon propagator 

We convert the ellipsoidal bound on configurations, 
just obtained, to a bound on the gluon propagator by 
taking expectation values, 



r 6c /-Vv p^pu E 



< 



k;\k+p\>\q\ 



N 2 q 2 



q 2 — p 2 



(k + p) 2 — p 2 
+ (2t:/L)J2 



(Bl) 
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From the fouricr expansions, 

(4(^)^(0)) = V-^S^D^We*" 

k 

= ( B2 ) 



where D^ik) is the gluon propagator on the finite peri- 
odic volume V, we obtain 



(4>U = v-H hd D^{k), 

and the last bound becomes 

P^D /J , v (k)p u 



(B3) 



v- 1 £ 



k; \k+p\ >\q\ 



(k + p) 2 — p 2 



< 



Nq 2 



2 ? 

q z — jr 



(2tt/L)J2W\1 (B4) 



where we have used f abc f abe = N8 ce . This is an ex- 
act bound on the gluon propagator on a finite periodic 
Euclidean volume V = L d . 

We take the infinite- volume limit, L — » oo, keeping 
k, q, and p finite, and obtain 



v- 1 £ 



P^D^(k)p, y Nq 



k\ | k+p \ > | q 

We divide out a factor of p 2 



(k + p) 2 — p 2 q 2 — p 2 



(B5) 



v- 1 E 

k; \k-\-p\ > |q| 



VlJ,D^{k)p v Nq 2 



(k + p) 2 — p 2 q 2 — p 2 ' 



(B6) 



where p M = p^/\p\ is a unit Lorentz vector. Recall that 
|p|,j3"/i and q^ are quantities at our disposal. We take the 
limit \p\ — > 0, keeping p^ and q fixed, which gives 



y- 1 E 

k;\k\>\q 



PnDfi V {k)p v 
k 2 



< N. 



(B7) 



We now take the limit q — > 0, and convert the sum to 
an integral, since we are in the infinite- volume limit, and 
obtain the bound 



1 



(2n)° 



d d k SiggW < N . 



(B8) 



Spherical symmetry is regained in the infinite-volume 
limit so, by transversality, D^ v {k) = D{k){8^ u — 
k^kv/k 2 ), and the angular average over k yields the 
bound 



d-l 1 

d {2n) d 



k z 



< N. 



(B9) 



Appendix C: Simple ellipsoidal bound on continuum 
configurations at infinite volume 

We note parenthetically that we could have taken the 
limit of large volume V, without taking expectation val- 
ues. In (|A33|) . we take the limit L — >• oo keeping k, q and 
p finite, divide by p 2 , take the limit p —> 0, followed by 
the limit q — > 0, and we obtain the ellipsoidal bound 



f abc f ade rf*rfp^Y, 



b* d 
u k,fi u k,u 



k 2 



(CI) 



We obtain a simpler ellipsoidal bound by summing over 

a complete basis EfeLi" 1 Vb* 7 lt = S ce , and Y,\=iPjlPu = 
5^, which gives, 



E 



b* n b 



< N(N 2 - l)d. 



(C2) 



In position space this bound reads, by ©, 



d d x A b Jx)[{-df l )- 1 A] b Jx) < N(N 2 - l)dV. (C3) 



For a configuration A Ax) that has compact support, this 
bound becomes vacuous in the limit V — > oo bcause the 
right-hand side diverges with V, while the left-hand side 
remains finite. However for a typical gauge-fixed config- 
uration (and here we may suppose that a lattice regular- 
ization is in place) , the integral on the left is a bulk quan- 
tity of order V , and the bound is meaningful. Bounds on 
lattice configurations are given in Appendix D. 



Appendix D: Bound on lattice gluon propagator 

1. Notation for lattice quantities 

Lattice configurations are defined by link variables 
U Xilx £ SU(N), that live on the link (x,x + e M ), where 
sites of the lattice are labeled (in lattice units) by inte- 
ger x v , and e M is a unit Lorentz vector in the positive fx 
direction. Numerical gauge fixing is done by minimizing 
the function 



J^07)=5^Retr(l-^ 



(Dl) 



with respect to local gauge transformations g x , where 
9 U x>fi = g~ l U Xtfi g x + efl is the gauge transform of the con- 
figuration U XtfM by g x . In practice there are many local 
minima, and the particular minimum chosen is algorithm 
dependent. For our purposes, any minimum will do; the 
absolute minimum plays no special role. The only prop- 
erties we shall use are that at any local minimum (i) the 
functional Fjj{g) is stationary, and (ii) the matrix of its 
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second derivatives is positive. Property (i) gives the lat- 
tice transversality condition, 

E^-^- e ^)=0, (D2) 
ft 

where we have introduced the lattice variables, 

A b x jU)^-tv[t\U x ^-Ul^], (D3) 

that are a lattice analog of the continuum variables 
A b (x). Property (ii) is the positivity of the matrix el- 
ement 

(V>,M(C7)V)>0, (D4) 

for all Here M^(U) is the lattice Faddeev- Popov 
matrix. It is a real symmetric matrix that is conveniently 
expressed as 

M = Mq + Mi, (D5) 
where Mq and Mi are defined by the quadratic forms 

X.fl 

X(^ + ^)(^+e M -^) ], (D6) 

where ij) x = t a ip x and ijj* = t a ip%*> an< ^ 

(^M 1 (A)^) = tr{[^ +eM ,^] (U^-U^)} 

= -(i/2)r bc (^ +eM +^r 

xA^(^ +6fl -Vx) c .(D7) 

The relation to the continuum Faddeev-Popov opera- 
tor ([3]) is apparent. The last expression is real when A 
satisfies the lattice transversality condition (|D2|) . Prop- 
erties (i) and (ii) define the (first) lattice Gribov region O. 

We consider a hypercubic periodic lattice of volume 
V = JV d , where A/" is an integer and = 0, 1, ...TV — 
l,mod(AQ. Fourier transformation is given by 

<^E a L cx P[*M^ + ^/2)] (D8) 

k 

where fc M = and = 0, 1...W — l,mod(7V). 

The transversality condition (|D2j) is diagonal in momen- 
tum space where it reads 

E A >t=C (D9) 

and we have introduced 

K„ =2sin(fc M /2), (D10) 
and similarly P M = 2sin(p M /2), = 2 sin(<7 A1 /2). 



2. Ellipsoidal bound on lattice configurations 

We proceed as in the continuum case. The lattice Gri- 
bov region Q is defined by the condition on (transverse) 
configurations U , 

- (V>,Mi(A)V>) < &,M (U)i>), (Dll) 

for all 0, where A = A(£^)- 

As in the continuum case, the matrix Mi (A) is linear in 
A. However in the lattice case, Mq{U) is not independent 
of the configuration U. We nevertheless obtain a simple 
lattice bound by introducing the matrix /C defined by 

{il),K i>) =^tr[ - 2{i> x+eti - ip x )* {ip x+e „ -ip x ) }, 

(D12) 

which is independent of U [22[. The difference, 

[/c - M (um = y, tr[ - «s+ e „ - r x ) 

X,fJ, 

x(l - U x ^)(l - ^V)(^ +e , - V.) ], (D13) 

is manifestly positive for every lattice configuration U 
and every trial wave function tp, so we have 

{i),M {A)i>)<{i>,1C i>) (D14) 

which, by (|Dllj) . yields the inequality 

-{^M 1 {A)i,)<{i>,K^), (D15) 

for every configuration U € £1 and every ip. Geometri- 
cally, it is natural to define a region O in configuration 
space by this condition, 

9 = {U : -{ip, M 1 {A)%1>) < (tp, K-o^) for all 0}, (D16) 

where A = A(U) is transverse, and we have the inclusions 

A C O C 6. (D17) 

Here A is the fundamental modular region which consists 
of every configuration that is the absolute minimum of 
the minimizing function on its gauge orbit. We shall de- 
rive a bound on the lattice gluon propagator that holds 
for all transverse configurations in O, which then holds 
a fortiori for all configurations in the Gribov region fl. 
Because /Co is independent of U and because Mi (A) is 
linear in A = A(U), the proof goes just as in the contin- 
uum case, but with Mo — > /Co. 
We define 

0o = V- 1/2 e lp - x r), (D18) 

where V = Af d and = is a lattice momentum, 

and we have 

/C o 0o = P 2 4>o, (D19) 
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where P M is defined as in (|D10[) . The continuum proof 
goes through, with the substitutions p 2 —> P 2 , k 2 — > 
K 2 ,q 2 — > Q 2 , and M — > /C , and (|A4|) becomes 

K 2 >Q 2 >P 2 > 4 sin 2 (71- /AT). (D20) 

By (|D7|) . equation (|A31[) gets replaced by 

(D21) 



(IA26j) gets replaced by 

X - P 2 - » £ P M co S (^/2)rV*aLo,^ (D22) 



(1X271) by 



fabc a* b c 



< 2tan(7r/A/"), 



(D23) 



(1X291) by 



X<P 2 + 2tan(7r/A/')^]|sinp Al |, (D24) 



(M) by 



F < 



iV 2 Q 2 



Q 2 -P 2 
and finally (|A"33|) by 



P 2 + 2tan(7r/J\T) V] | sinp M 



(D25) 



k 

< N 2 Q 2 



Q 2 -P 2 



K' 2 - P 2 
\p 2 + 2 tan(?r/A0 ' sinp ^ 



(D26) 



where W„ = (fc+p)„, K' 2 = 4 £ M sin 2 (fc'/2), C7 M (fc',p) = 
P^ cos(/c^/2), and the sum over k is restricted by if' 2 > 
Q 2 . This is an ellipsoidal bound that holds on a finite 
lattice for every configuration a b k in the lattice Gribov 
region tt. 



3. Bound on lattice gluon propagator 

Upon taking expectation values, we obtain a bound on 
the lattice gluon propagator, 



k 



C^k',p)D^(k)C l/ (k',p) 
K' 2 - P 2 

NQ 2 



(D27) 



< 



- P 2 



P 2 + 2 tan(7r/A0 ^ | sinp M | 



where D^ u (k) is the gluon propagator on a finite lattice, 
and we have used 

« 4,u) = V- l 5 bd D,„(k). (D28) 
We now take the lattice volume to infinity, Af — > 00, 
while keeping the lattice spacing finite. The lattice 
momentum = l-Kn^jM becomes a continuous angle 
— 7r < fc M < 7T. We divide out P 2 , take the limit p — > 
and then q — > 0, and we obtain the bound on the gluon 
propagator on an infinite lattice, 



PuPu.P, < N, 



where P is an arbitrary unit vector. Here 



(D29) 



T 



(2ir) d 



d cos(A; M /2) D^(k)cos{k v /2) 



4£ A sin 2 (fc A /2) 



(D30) 



(no sum on fi or ^) is a tensor that is invariant under 
the hypercubic symmetries. It is thus of the form T^ v = 
J&iiv^ where J = T^^/d, and the bound on the lattice 
gluon propagator reads J < N, or 



J £, cos 2 (fc il /2) P>uu(fc) 
d d fc ^ M „ 1 M/ „,, , MMV 7 < AT. (D31) 



d(27r) c 



4£ A sin 2 (fc A /2) 
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